High cooperativity in coupled microwave resonator ferrimagnetic insulator hybrids 
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We report the observation of strong coupling between the exchange-coupled spins in gallium-doped 
yttrium iron garnet and a superconducting coplanar microwave resonator made from Nb. The mea- 
sured coupling rate of 450 MHz is proportional to the square-root of the number of exchange-coupled 
spins and well exceeds the loss rate of 50 MHz of the spin system. This demonstrates that exchange 
coupled systems are suitable for cavity quantum electrodynamics experiments, while allowing high 
integration densities due to their extraordinary high spin densities. Our results furthermore show, 
that experiments with multiple exchange-coupled spin systems interacting via a single resonator are 
within reach. 
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The study of the interaction of matter and light on the 
quantum level is at the core of solid state quantum infor- 
mation systems. Strong [HH] and ultra-strong coupling 
[2l[4] has been achieved, allowing for the coherent transfer 
of quantum information. For the practical implementa- 
tions of quantum information systems, the use of hybrid 
systems has been suggested. In such hybrids, natural mi- 
croscopic systems (atoms, molecules, electron spins, and 
nuclear spins) are coupled with artificial meso-scale struc- 
tures such as superconducting quantum circuits by means 
of microwave photons. Whereas the former have long 
coherence times due to sufficient decoupling from envi- 
ronmental noise, the latter allow for fast qubit gates due 
to strong coupling to electromagnetic fields [5]. While a 
large variety of hybrid systems has been proposed [SHE] , 
ensembles of electron spins as quantum memories [9] [10] 
seem most promising and their coupling to supercon- 
ducting resonators [TlTfTT] and flux qubits [TB] has been 
studied recently. Although the coupling strength g of 
an individual spin to an electromagnetic mode of a su- 
perconducting microwave resonator is small (typically in 
the order of 10 Hz), the coupling of an ensemble of N 
spins is enhanced by a factor of \/N \T§\ [20]. In this 
way, strong coupling g c s = gyN 3> /c, 7 can be realized, 
where k and 7 are the loss rates of the resonator and 
spin system, respectively. With loss rates in the order of 
MHz, typically 10 12 spins are needed to reach the strong 
coupling regime. Therefore, the spin system has to fill a 
major part of the small mode volume of a superconduct- 
ing resonator. Until today, mostly paramagnetic systems 
consisting of ensembles of noninteracting spins have been 
studied. The coherent coupling of microwave resonators 
to ferromagnetic systems with strongly exchange coupled 
spins remains to be explored. Soykal and Flatte [2TJ HI] 
theoretically discussed the strong coupling of photonic 
and magnetic modes in exchange locked ferromagnetic 



systems. One particular advantage of, ferromagnetic sys- 
tems is their higher spin density. Therefore, for the same 
number N of spins, their volume can be reduced consid- 
erable compared to paramagnetic systems. This should 
allow one to couple multiple spin ensembles to the same 
microwave resonator. 

In this letter, we investigate the coupling between the 
electromagnetic modes of a superconducting coplanar 
waveguide microwave resonator and the magnetic modes 
of the exchange-locked ferrimagnet yttrium iron garnet 
(Y 3 Fe 5 0i2 or YIG) doped with gallium (YIG:Ga). We 
measure a coupling rate of <? c ff/27r = 450 MHz exceed- 
ing both the spin relaxation rate 7/2-77 — 50 MHz and 
the resonator decay rate k/2it = 3 MHz. That is, we 
observe strong coupling. The measured effective cou- 
pling strength follows g e g = gy/N, where the number of 
spins interacting with the resonator is estimated from the 
sample geometry. Furthermore, the measured relaxation 
rate of the spin system is fully consistent with the nat- 
ural linewidth of YIG:Ga obtained from ferromagnetic 
resonance (FMR) measurements. We note that YIG is 
one of the prime candidates for studying strong coupling 
between exchange locked spins and the electromagnetic 
modes of a microwave resonator. It is known to have ex- 
cellent microwave properties, in particular a very small 
FMR linewidth of 10 pT at 4K and uj/2it = 9.3 GHz 
|23j . This narrow linewidth corresponds to a T2 time in 
the order of microseconds [23]. Since high quality YIG 
thin films can be prepared on various substrates (gadolin- 
ium gallium garnet [55J[55], Si and GaAs [57]) and doped 
with rare earth elements in order to adjust the FMR 
linewidth in a controlled way. YIG appears like an ideal 
material for ferromagnet based quantum hybrids. 

As pointed out by Soykal and Flatte [21] [22], in 
a macrospin approximation the Hamiltonian for the 
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FIG. 1. Schematic of the experimental setup. The (purple) 
gallium doped YIG sample is cemented ontop of one of the 
niobium microwave resonators which are arranged to allow 
for multiplexing. Experiments are performed at millikelvin 
temperatures in transmission by vector network analysis in a 
superconducting solenoid magnet. 



ferromagnet-resonator system can be expressed as 

H — Tiu) r a) a + — ~S Z + hg(aS+ + ($S-). (1) 

Here, a' and a are the photon creation/annihilation 
operators, uj r the resonator frequency, g s is the elec- 
tron g-factor, [1b the Bohr magneton, and B z the ap- 
plied magnetic field establishing the resonance condi- 
tion Uujfmr = 9s^B z /h between the microwave field 
and the energy level splitting of the ferromagnet. The 
macrospin operator S = (S+e- — S'_e + )/v / 2 + S z z with 
e± = =f(x ± «y)/v2 is expressed in terms of the spin 
lowering and raising operators 
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Here, \£, m) are the eigenstates of the macrospin and we 
have assumed that the macrospin state is fixed at its 
maximal value £ = N/2. In the Dicke model [25J of N in- 
dependent paramagnetic spins this would correspond to 
a fully excited spin system with no photons in the cavity. 
In contrast to the Dicke model, for our macrospin model 
states with £ < N/2 are not accessible. Due to strong 
exchange coupling, states with £ < N/2 are separated 
in energy and require the excitation of magnons. We 
note that the coupling between the photonic and mag- 
netic system is a magnetic dipole transition and that the 
Hamiltonian conserves the total excitation number 
Z = n + m, where n is the photon number in the cavity 



and \m\ < £ — N/2 the magnetic quantum number. As- 
suming that S is antiparallel to B z and n = initially, 
we have Z = N/2 and, hence, can index the basis states 
1 77, m) of the resonator-spin systems either by the photon 
number 77 (|n, ^ — n)) or the magnetic quantum num- 
ber to (|y — to,™)). Evidently, these basis states are 
similar to those of the Dicke model [28] for a paramag- 
netic ensemble of N noninteracting spins coupled to a 
resonator, with £ = N/2 taking the role of the cooper- 
ation number. Note that the coupling between the fer- 
romagnet and the resonator modes changes over a wide 
range when moving through all possible states \n, -g- — 77) 
with an equilibrium point at 770 = (2N — l)/3, where the 
ferromagnet-resonator systems shows superradiance [35] . 

Due to the analogy with the Dicke model, we expect 
that the coupling strength of the ferromagnet-resonator 
system is given by the effective coupling strength g e s = 
gy/N of a paramagnet-resonator system, where g = 
Bifi is the coupling rate of an individual spin with 
the magnetic quantum number 777 = 1/2 to the resonator 
|28) . It is determined by the magnetic component of the 
rf vacuum field Bi q = yf pQhuj c /2V lrL which depends on 
the resonance frequency of the microwave resonator uj c 
and the mode volume V m of the resonator 29J . Hence, for 
a given resonance frequency, g c s = 9s ^ l B \J pophoj 1 V/2V lrL 
depends only on the spin density p = N/V and the filling 
factor V/V m , where V denotes the volume of the overlap 
between the spin system and the resonator field mode. 
Non-interacting spin ensembles like paramagnetic centers 
in semiconductors or insulators typically have a spin den- 
sity p in the order of 10 15 < p < 10 18 cm" 3 [TT 1 [H ) [T3] . 
In combination with resonator frequencies in the 5 GHz 
range, this results in coupling rates in the order of 10 MHz 
assuming V/V m ~ 1. In contrast, exchange coupled sys- 
tems naturally have a spin density in the order of one per 
atom (e.g. Fe, Ni, Co) or in the case of YIG 40 per unit 
cell (unit cell volume: a 3 = 1.8956 nm 3 ), corresponding 
to a spin density of 2 x 10 22 cm -3 [30] ■ Due to the increase 
of at least four orders of magnitude in the spin density we 
expect in exchange coupled systems a two orders larger 
coupling strength as compared to non-interacting spins 
using similar microwave resonators. In other words, in 
exchange coupled systems the sample volume of the spin 
system can be reduced by a factor of 10 4 while keeping 
the coupling rate constant, thereby allowing for a much 
higher integration density. 

In our experiments, we study the coupling between the 
exchange locked system YIG:Ga and a superconducting 
Nb resonator. The resonator structure is patterned into 
a 100 nm thick Nb film deposited onto an intrinsic sili- 
con substrate using optical lithography and reactive ion 
etching [29] . Figure [l] shows the layout of the microwave 
circuitry consisting of an input line, three resonators with 
resonance frequencies of Ja = 5.90 GHz, f B = 5.53 GHz, 
and fc = 5.30 GHz, and an output line. This configura- 
tion allows to compare loaded and unloaded microwave 
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FIG. 2. Transmission spectrum of the YIG-microwave res- 
onator hybrid taken at T = 50 mK as a function of the applied 
magnetic field B z . The coplanar waveguide resonators (A-C) 
show a slightly decreasing resonance frequency with increas- 
ing in-plane magnetic field. Additionally, resonator A shows 
a pronounced avoided crossing at B z ,r = 170 mT, where the 
resonance condition Tiu r = g 3 fisB ZiI /h between the resonator 
and the energy level splitting of the ferromagnet is satisfied. 
The red line shows a fit according to eq.(|3|. 



resonators on the same chip. A 2 x 0.5 x 0.7mm 3 sized 
(length x width x thickness) commercial YIG:Ga crystal 
is cemented onto resonator A with the highest microwave 
frequency j a = 5.90 GHz. The number of spins interact- 
ing with the resonator is estimated from the overlap of 
the YIG:Ga crystal with the meandering coplanar waveg- 
uide with a center conductor width of 6 u.m and a gap of 
12 u.m. With the overlap length of 2.5 mm and assuming 
that the vertical extension of the microwave field into the 
YIG crystal is about 30 u.m, the total number of spins 
coupled to the resonator is estimated to N ps 4.5 x 10 16 . 
To preserve the superconducting state of the microwave 
resonators, the surface of the chip is carefully aligned in 
parallel to the applied magnetic field B z generated by a 
superconducting solenoid. The microwave transmission 
experiments are performed at the base temperature of a 
dilution refrigerator of 50 mK using a commercial vector 
network analyzer. To thermalize the microwave input 
and output lines, attenuators are used at the 4 K, the 
still and the mixing chamber stages. 

Figure [2] shows the microwave transmission IS21I 2 as a 
function of frequency and applied magnetic field. Hereby, 
52i is the complex scattering parameter. In the spec- 
trum at B z = 0, four transmission peaks are visible cor- 
responding to the resonance frequencies of the coplanar 
microwave resonators A, B, and C. The broad feature 
labeled D stems from a parasitic mode of the metallic 
microwave box in which the sample is mounted. As ex- 
pected, the resonators B and C show only a weak mag- 
netic field dependence, because they are not interacting 
with the YIG:Ga crystal due to the absence of physical 



overlap (cf. Figs. [I] and [2]). On the contrary, resonator 
A and the box mode D couple to the magnetic system. 
While mode D allows us to probe the ferromagnetic res- 
onance independently of the strongly coupled mode (A) 
and to determine the resonance condition, resonator A 
shows a distinct anticrossing at B ZT — 170 mT where 
the resonance condition Hco T = g s ^BB z ^/h between the 
resonator frequency and the energy level splitting of the 
ferromagnet is satisfied. 

To derive the effective coupling rate g c g from the mea- 
sured data, we model the ferrimagnet-resonator interac- 
tion according to two-level atom theory. Within this ap- 
proach, the dispersion of the resonance frequency is given 

by m 



(3) 



Here, A = w r — wfmr — 9s^B{B ZtT — -Bfmr)/?* is the field 
dependent detuning between the resonator frequency ui r 
and the field dependent FMR frequency wfmr- The ex- 
perimental data agree very well with this model pre- 
diction. Fitting the data yields g ff = 450 MHz and 
Bfmr — 170 mT, using the (/-factor of the ferrimagnetic 
resonance g s — 2.17. For the geometry of our resonator 
we estimate g/27r ~ 5 Hz and N ~ 4 x 10 16 . With these 
numbers we expect g c g/2n = (g/2n) VN ~ 1 GHz cor- 
roborating our experimental result within a factor of two. 
We point out that no gap has been assumed between 
the microwave resonator and the YIG:Ga sample in our 
estimate of g and N. This most likely is not the case 
in the actual sample configuration due to imperfections 
in the sample mounting process. Here, an additional 
gap between the superconducting niobium layer and the 
YIG:Ga would result in a decrease in N. We finally note 
that the frequency dispersion of the ferromagnetic mode 
compares well with values from Ref. [32], whereas the g- 
factor of 2.17 is about 10 % larger than g s = 2.0 reported 
for room temperature YIG [35] . This is attributed to an 
imperfect calibration of the magnetic field produced by 
the superconducting solenoid. 

To get information on the relaxation rate 7 of the spin 
system we analyze the evolution of the linewidth of the 
resonator mode A as a function of the magnetic field 
B z by fitting a Lorentzian lineshape in the frequency 
domain for every measured magnetic field magnitude 
[17) . Figure |3][a) shows the resonance frequency obtained 
from such a fit as red crosses superimposed on the color- 
coded dataset. In Fig. [3^b) the corresponding (FWHM) 
linewidth (red crosses) is shown. At low magnetic fields, 
the resonator mode A is essentially decoupled from the 
spin system. In this regime the measured linewidth is 
given by the resonator loss rate k. Closer to the ferro- 
magnetic resonance, the linewidth of the system is given 
by the combined relaxation rate of the spin system and 
the microwave resonator leading to an increase in the 
observed linewidth. 
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FIG. 3. Analysis of the resonance frequency and linewidth as 
a function of the magnetic field B z . Panel (a)shows the res- 
onator transmission |5*2i| 2 (same data as Fig.[2J) as a function 
of frequency and applied magnetic field. The red crosses mark 
the resonance frequency determined by fitting a Lorentzian 
lineshape to the data at constant B z . The red crosses in panel 

(b) show the extracted FWHM linewidths. In addition, the 
blue circles show the linewidths obtained from the numeri- 
cal simulation of the transmission spectra plotted in panel 

(c) . The simulation is based on the input-output formalism 
resulting in eq.(g[) pH H7| [20l |33] . 

We now analyze the measured data in more detail to 
quantify the coupling and loss rates in our system. Here, 
we use a standard input-output formalism [H | ll7 1 l2"C)] , l33|. 
Within this framework the transmission of microwaves 
from the input to the output port of the microwave res- 
onator is given by 
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Here, u>/2n is the frequency of the microwave probe tone, 
k c is the external coupling rate between the microwave 
resonator and the feed line, and k; summarizes the intrin- 
sic loss rate of the microwave resonator. In our case we 
have Ki <C k c , resulting in a total microwave resonator 
relaxation rate n ~ k c . Fig.[3]jc) shows the calculated 
transmission using g c g/2ir — 450 MHz, j/2ir — 50 MHz, 
and k/2tt — 3 MHz. Evidently all features of the exper- 
imental data of Fig.[3]ja) are nicely reproduced. In par- 
ticular, Fig.[3]jc) clearly shows that for the large cavity 
decay rate k/2tt — 3 MHz the two transmission peaks ex- 
pected at B rcs = 170 mT cannot be resolved due to the 
limited signal to noise ratio in the experimental data. 
Moreover, we can analyze the simulation data shown in 
Fig. 3 c) in the same way as the experimental data in 
Fig. 3 a) to derive the field dependence of the linewidth. 
The result is shown by the blue circles in Fig.|3jb), where 
in contrast to the experimental data, the modeled data 
is noise-free allowing to predict the linewidth for all mag- 
netic field values. The good agreement between experi- 
mental and simulation data again demonstrates that the 



parameters chosen in the simulation well reproduce the 
experimental situation. In summary, our analysis shows 
that g c g ^ k, 7 with a cooperativity C = g^ s / K7 ~ 1350. 
That is, the strong coupling regime has been reached for 
the ferrimagnet-resonator system. 

Finally, we compare the experimentally determined re- 
laxation rate 7 with measurements of the temperature de- 
pendence of the FMR linewidth performed at 9.43 GHz. 
Rachford et al. jM] report linewidth values for lightly Ga 
doped YIG that decrease from 1 mT at liquid helium tem- 
peratures to 0.1 mT at room temperature, corresponding 
to 28 MHz and 2.8 MHz, respectively. Our sample has a 
higher Ga-doping concentration [35] . which is known to 
result in an increase of the linewidth. This is consistent 
with the relaxation rate of 7/27r = 50 MHz observed in 
our experiments at millikelvin temperatures. Although, 
we would like to stress that little is known about the 
evolution of damping in YIG for T < 2K, according to 
Sparks and Kittel [36] the limiting relaxation mechanism 
is spin-lattice coupling, which is expected to be well be- 
low 1 MHz for YIG. 

In conclusion, we experimentally observed strong cou- 
pling between a superconducting microwave resonator 
and a gallium doped yttrium iron garnet ferrimagnet. 
The effective coupling rate of 450 MHz reaches 13% of 
the resonator frequency w r and by far exceeds the relax- 
ation rate 7/27r = 50 MHz of the spin system. Therefore, 
exchange coupled spin systems indeed can be used for 
cavity quantum electrodynamics. Furthermore, the large 
coupling rates achievable in exchange coupled systems al- 
low to place more than a single magnetic system in the 
microwave resonator. This should allow to study e.g. the 
exchange of magnetic excitations via a cavity bus similar 
to the approaches pursued in the field of cavity QED [5] . 
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